Transport Properties of a Delta-Shell Gas with Long Scattering Lengths 
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The coefficients of diffusion, thermal conductivity, and shear viscosity are calculated for a system 
of non-relativistic particles interacting via a delta-shell potential V{r) = —v6{r — a) when the 
average distance between particles is smaller than a. The roles of resonances and long scattering 
lengths including the unitary limit are examined. Results for ratios of diffusion to viscosity and 
viscosity to entropy density are presented for varying scattering lengths. 
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Gottfried's treatment [T, ? of the quantum mechanics 
of two particles interacting through a delta-shell poten- 
tial V(r) = —vS{r — a), where v is the strength and a is 
the range, forms the basis of this work in which the trans- 
port properties of a dilute gas of non-relativistic particles 
are calculated for densities n for which the average inter- 
particle distance d ~ n~^/^ a. The delta-shell po- 
tential is particularly intriguing as the scattering length 
asi and the effective range tq can be tuned at will, in- 
cluding the interesting case of asi oo that is accessible 
in current atomic physics experiments (see, for e.g., [3 ). 
The delta-shell potential also allows us to understand the 
transport characteristics of nuclear systems in which the 
neutron-proton and neutron-neutron scattering lengths 
(asi — — 23.8 fm and —18.5 fm, respectively [J) are much 
longer than the few fm ranges of strong interactions. In 
addition, the delta-shell potential admits resonances. 

For small departures of the distribution function from 
its equilibrium, the net flows of mass, energy and mo- 
mentum are characterized by the coefficients of diffusion, 
thermal conductivity and shear viscosity, respectively. 
Here we calculate these coefficients in first and second or- 
der of deviations using the Chapman-Enskog method [5 . 
In particular, we examine these transport properties for 
varying physical situations including the case of an infi- 
nite scattering length, i.e., the unitary limit, and when 
resonaces are present. 

We begin by collecting the two-particle quantum me- 
chanical inputs [ij |2] [6] and cast them in forms suit- 
able for the calculation of the transport coefficients. The 
Hamiltonian for two particles interacting via the delta- 
shell potential is = —li?^—v5{r — a) , where ji denotes 

the reduced mass, A = \^r'^^ — ^ is the Laplacian 

" or or ^ 

in spherical coordinates (L is the orbital angular momen- 
tum operator), r is the separation distance, v and a are 
the strength and range parameters of the potential, re- 
spectively. In terms of the dimensionless variable p = kr 
{k is the wave number), E = and / the orbital quan- 
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tum number, the general solution of Shrodinger equation 
is given by the spherical Bessel functions ji{p) and ni{p): 

^{p) = u{p)/p = Aijiip) + Bimip) • (1) 

The phase shifts 5i are obtained through [T] 

gxjf{x) 



tan(J^ 



(2) 



1^ gxji{x)ni{x) ' 

where x = ka^ and the single dimensionless parameter 

g = {2fiv/h^) a (3) 

controls all physical outcomes. At low energies (/ = 0), 
the scattering length a^/, the effective range ro and the 
shape parameter P, are defined through |7j : 

kcot{So) = - — + ro^ -Prlk^ + 0(A:^) , 

dsl ^ 
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Note that as ^ ^ 1, a^^ ^ oo, but ro/a ^ 4/3 and P 
3/80. For / > 0, Newton's generalization of scattering 
lengths and effective range parameters reads as |7] 



cot{Si) 



1 



.(0 



(0 +r^^Y+0(fc^). 



For the delta-shell potential, the I > results are 
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The scattering lengths in Eq. (p) become very large for 
values of ^ ^ 2/ + 1. In Fig. T} we show asi and ro 
as functions of the strength parameter ^ for / = 0. The 
diverging tendencies of of |as/| as ^ ^ 1 and |ro| as ^ ^ 
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FIG. 1: (Color on-line). Scattering lengths and effective 
ranges as functions of the strength parameter g in Eq. (|3|. 



are clearly evident. The case of ^ — oo corresponds to 
the hard-sphere case. 

We turn now to calculate the coefficients of diffusion, 
thermal conductivity and shear viscosity, based on the 
Chapman-Enskog approach [5 , 8 . The central ingredient 
is the transport cross-section of order n: 



j dcos6>(l -cos^6>) 



dft I 



(7) 



where the scattering angle and the collisional differen- 
tial cross-section ^^IMll ^^^e calculated in the cen- 

cL^l \c.m. 

ter of mass reference frame of the two colliding particles. 
For indistinguishable particles, an expansion of the cross- 
section in partial waves Xl^(2/ + l)(e*^'^^ — l)P/(cos^) and 
the orthogonality of the Legendre polynomials Pi renders 
the integral above to the infinite sums 



,(1) 



47ra2 



,(2) 



47ra2 



X 



{21 + 3) 



^'(2? + l)sin2(,50, 

1+2 -Si) (8) 



The prime on the summation sign indicates the use of 
even / for Bosons and odd / for Fermions. In Eqs. ([8|, the 
low energy hard-sphere cross section {g — oo) 47ra^ has 
been used to render the transport cross sections dimen- 
sionless. If the particles possess spin 5, then the properly 
symmetrized forms are: 



QFermi^ for integer s , 



(n) _ g + 1 (n) 

= ^ 2^Tl ^""'^ half-integer s . 

Here, we will present results for the case of spin-1/2 par- 
ticles only. The transport coefficients are given in terms 
of the transport integrals 



cj^^'*^(T) = / ^7 6-"^ Y^^q^^'^x), (9) 
Jo 



where 7 = 



hk 



( KT) 



de-Broglie wavelength A = (27rfi^/m /c^T)^/^. 



with the thermal 



In what follows, the coefficients of self diffusion ^, 
shear viscosity 7^, and thermal conductivity tz are nor- 
malized to the corresponding hard-sphere-like values 



3V2 h 

32 mna^ 



T] 



5V2 h 

"32" 



and 



75 hkB 
64^2^ 



(10) 

In the first order of deviations from the equilibrium dis- 
tribution function, the transport coefficients are 

1 



Ml 

V 



A(T);^(M)(^)' 

1 



Equation (11) shows clearly that if 



,(2,2) 



(11) 



is T- 



independent(as for hard-spheres with a constant cross 
section), the shear viscosity exhibits a T^/^ dependence 
which arises solely from its inverse dependence with A(T). 
For energy-dependent cross sections, however, the tem- 
perature dependence of the viscosity is sensitive also to 
the temperature dependence of the omega-integral. 
The second order results can be cast as 



= (l^5^(T)){l±n\\^(T)) , 



(12) 



where is ^ or 7^ or and the ± refers to Bose (+) and 
Fermi (— ) statistics. For example. 



(2,2) 



2uj- 



(2,3)x2 



ix) 



-7/2 



(2,2) 



(2,2) 
128 ^4/3 

33/2 J2,2) 



(13) 



It is useful to define a characteristic temperature 



T 



or 



(14) 



kBTna? T 

in terms of which limiting forms of the transport coef- 
ficients can be studied. For spin 1/2 statistics. Fig. 2^ 
shows from Eq. ( 11 ) normalized to & from Eq. (^10). 
As expected, the diffusion coefficient grows steadily with 
temperature. For g ^ 1,3, and for T <t^T^ we find the 
asymptotic behavior (the dashed lines in Fig. [2| 

2( 1 a/I- (15) 
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The cases ^^1,3 require special consideration. In these 
cases, the asymptotic behavior for T <C T is obtained 
from a series expansion in x = ka of the transport cross 
sections q^'^\x) in Eqs. ([8|, and subsequent integrations 
of Eq. ([9|. We find the limiting forms: 



Stt 
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for g 
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FIG. 2: Normalized (to ^ in Eq. (10)) diffusion coefficient as 
a function of normalized (to T in Eq. (14)) temperature for 
various values of the strength parameter g in Eq. Solid 
curves show the diffusion coefficient in the first approxima- 
tion, Vi in Eq. (11), and the dashed lines are its asymptotic 



trends for T < T in Eq. ( 14) 
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FIG. 3: (Co lor on-line). The normalized viscosity coefficient 
in Eq. as a function of the inverse scattering length. 

Effects due to resonances associated with the partial waves 
/ = and 1 are indicated by the letters S and P, respectively. 
In the vertical shaded region, a large number of partial waves 
are required to obtain convergent results. 



The asymptotic behavior (for T <C T) of the shear vis- 
cosity [t]]-^ is similar to that of the diffusion coefficient: 



7/ 
fj 



(¥) (^/^) for ^^1,3 

/ ~\3/2 

67r \T/T\ ioi g = l 

^ 3^(T/f)'^' for ^ = 3. 



(17) 



It is interesting that even at the two-body level, the co- 
efficients of diffusion, thermal conductivity and viscosity 
acquire a significantly larger temperature dependence as 
the scattering length a^^ ^ oo (^ ^ 1). 

Figure [3] shows the normalized viscosity coefficient as 
a function of — = for T < T. Enhanced cross 

asi g 

sections at resonances produce significant drops in the 
viscosity as g 2/ + 1. The widths of the dips in viscosity 



decrease with increasing values of /. For T < the dips 
become less prominent and disappear for T >T. For the 
densities and temperatures considered, the contribution 
from the second approximation is small as < 0.05. 
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FIG. 4: (Color on-line). Ratio of diffusion (times mn) to 
shear viscosity versus inverse scattering length. In the vertical 
shaded region, a large number of partial waves are required to 
obtain convergent results. The hard-sphere results {g —00) 
are shown by solid circles. 
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TABLE I: Firsr order coefficients of diffusion (times mn), 
shear viscosity, and their ratios for T <^T for select strength 
parameters g. The unitary limit {g = 1) result for 77 was 
obtained earlier in Ref. 

In Fig.[4j the ratio of the coefficients of diffusion (times 
mn) and viscosity are shown as functions of a/asi for 
T < T. As expected, the largest variations in this ratio 
occur as ^ ^ 2/ + 1, that is, as resonances are approached. 
For T ^ T, the ratio approaches the asymptotic value 
6/5 away from resonances and 4/5 for the 5- wave res- 
onance (a point). With increasing T, resonances be- 
come progressively broader with diminishing strengths; 
for T ^ T resonances disappear. For comparison, this 
figure also includes results for hard spheres. 

The coefficent of viscosity (as also mn times the coef- 
ficient of diffusion) has the dimension of action (h) per 
unit volume. The manner in which the effective physical 
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volume V changes as the strength parameter g is var- 
ied is ihuminating as our results for T <C T in Table I 
shows. In the unitary limit = 1), the relevant volume 
is V (X A^. For ^ = 3, V cx Xa^ (independent of a^^), and 
for 1,3, V oc Xa'li (independent of a). 
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FIG. 5: (Color on-line). The insert shows the occurence of a 
minimum in the ratio viscosity to entropy density versus T. 
For each na^, this minimum value in the first (lower curve) 
and second (upper curve) order calculations of shear viscosity 
is plotted versus a/asi for the indicated na^ . In the verti- 
cal shaded region, a large number partial waves are needed. 
Vertical lines with letters {S,P,D, and F, respectively) indi- 
cate resonances associated with the partial waves / = 0, 1,2 
and 3. The symbol HS denotes hard-spheres result for which 
g ^ —oo. The horizontal dashed line shows the conjectured 
lower bound 1 / (47r) ^Oj . 

Recently, a lower limit to the ratio of shear viscos- 
ity to entropy density is being sought [11] with results 
even lower than {4:7r)~^ {h/kB) first proposed in [10]. We 
therefore examine 1,2/5 using the entropy density 

s = {5/2 -\n{nX^) ^ Ss{T)na^)nkB , 
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which includes the second virial correction [8^ to the ideal 
gas entropy density. The second virial coefficient [8] 

a2(T)=T2-^/^-23/2^'(2/ + l) 




where the prime indicates summation over even /'s for 
Bosons (— ) and odd Ts for Fermions (-h), Ei is the en- 
ergy of the bound state with angular momentum / and 
^(T) = (A/a)V(27r). The insert in Fig. [s] shows a char- 
acteristic minimum of rj/s versus T for fixed dilution na^ 
and strength g. Values of [77] 1, 2/^5 at the minimum are 
also shown functions of a/agi for several na^. The lower 
(upper) curve for each na^ corresponds to the first (sec- 
ond) order calculations of r]. The large role of improved 
estimates of r] on the ratio 77/5 is noticeable. The case 
of na^ = 1 possibly requires an adequate treatment of 
many-body effects not considered here. We can, how- 
ever, conclude that in the dilute gas limit r]/s for the 
delta-shell gas remains above {l/4:7T)h/kB- 

Our analysis here of the transport coefficients of par- 
ticles subject to a delta-shell potential has been de- 
voted to the dilute gas (non-degenerate) limit, in which 
two-particle interactions dominate, but with scattering 
lengths that can take various values including infinity. 
Even at the two-body level considered, a rich structure 
in the temperature dependence and the effective physical 
volume responsible for the overall behavior of the trans- 
port coefficients are evident. The role of resonances in 
reducing the transport coefficients are amply delineated. 
Matching our results to those of intermediate and ex- 
treme degeneracies [9^, T2^ which highlight the additional 
roles of superfiuidty and superconductivity reveals the 
extent to which many-body effects play a crucial role. 
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